Optical cooling and trapping of tripod-type atoms with rectified radiation forces. 
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A new scheme of three-dimensional (3D) all-optical (nonmagnetic) cooling and trapping of reso- 
nant atoms, based on using of so-called rectified radiation forces in non-monochromatic light fields 
is presented. It can be applied to the atoms with a tripod-type configuration of levels: atoms (ions) 
with the quantum transition F — 1 — > F = 0. 

The scheme proposed provides a long-term trapping of such atoms in deep light-induced potential 
wells. Moreover, the atom temperature can continuously be changed by varying field parameters in 
quite a large range (from super-Doppler to sub-Doppler values) without violating the localization 
stability. 
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I. INTRODUCTION 

The problem of optical cooling and trapping of atoms 
plays one of the central roles in studying resonant light 
pressure. [IHl] . The most famous useful solutions of this 
problem are the 3D confinement and cooling of atoms 
in a magneto-optical trap (MOT) [4[ and optical mo- 
lasses (OM) An indispensable element of MOT is 
a nonuniform magnetic field allowing one to "circum- 
vent" the optical Earnshaw theorem (OET) Q. In OM 
the principle of the 3D all-optical (nonmagnetic) viscous 
confinement of particles is manifested and, as contrasted 
to MOT, any confining space-dependent restoring force 
towards the center of OM is absent. In OM the vis- 
cous damping force acts upon the atoms and the damp- 
ing both of chaotic (resulting from laser cooling) and di- 
rected (macroscopic) atom movement occurs. Therefore, 
the time of diffusion escape of the atoms from the OM 
area may be quite long as compared to the time of free 
expansion of the atom cloud. Note that OM is an exam- 
ple of an all-optical (nonmagnetic) device, allowing one 
to perform 3D cooling of the atoms to ultralow temper- 
atures - below the so-called Doppler limit [3|. 

The number of efficient all-optical methods of cooling 
and trapping of atoms can significantly be increased when 
using rectified radiation forces (RRFs), induced by non- 
monochromatic (particularly, biharmonic ) nonuniform 
fields such as standing waves. RRFs were discovered and 
described in [7|, \m, and experimentally demonstrated for 
the first time in Q. Later, the idea of RRF was devel- 
oped in many other studies, see, for example, [l(| and 
the references given there. 

RRFs appear as the consequence of the nonlinear in- 
terference phenomena exhibited under the interaction of 
an atom with polychromatic electromagnetic field [8] . To 
explain their appearance physical mechanism we consider 
a simple model of the gradient force rectification, based 



on the results of the original works [7, 8] . Let a two- level 
atom moves with a velocity v in a strong bichromatic 
standing wave field. Local Rabi frequency has the form 



(V (x) + V 1 (x)e 



where Vq(x) = Vocos(kx), V\(x) = Vi cos[(fc + 6k)x + 
<p], Aq, Ai + Ao are the detunings from resonance and 
we assume the following hierarchy of the characteristic 
frequencies (7 is the spontaneous relaxation rate). 
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Then, for weak saturation V^/AiAo, Vq 2 /Aq <C I and 
5k = Ai/c <C k, the mean dipole force F (acting upon 
the atom) can be written as (more details see in [7, 8]): 



dx 
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dx 



dx 



where III and II2 are the populations of the adiabatic 
(dressed) atom states, Q = II2 — IIi is the population 
difference, U g = Vq(x)/A + Vi(x)/Ai is potential of 
gradient force F g = —dll g /dx determined by Stark (light) 
shifts of the atomic energy levels. This expression has a 
clear physical meaning: the atom in adiabatic state |1) 
moves in the force field with the potential U g and atom in 
state |2) moves in the force field with the potential (—U g ). 
Spontaneous transitions cause the incoherent mixing of 
the adiabatic states, moreover the transition rates be- 
tween the adiabatic states depend on the saturation pa- 
rameter V^(x)/Al [1, 7, 8]. 

The result of this dependence is the spatial modula- 
tion (with the period ~ A/2 = n/k) of the population 
difference Q = Q(x). In the adiabatic limit: kv/j —> 0, 
Q and the force F can be presented as (see [1, 8]) 
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where the term Fi = Q\(x)dVi(x) / dx one can interpret 
as the interference (nonadditive) contribution to the ra- 
diation force F (Fi = if V x = or V = 0). Evidently, 
that the gradient force F g = — (dU g /dx) does disappear 
(after averaging over the spatial oscillations with micro- 
scopic periods ~ 1/fc), but the force Fi does not disap- 
pear (as the change in sign of the gradient of the Vi(x) 
can be compensated by the change in sign of the oscillat- 
ing component of the Q\ due to phase shift VP between 
the two standing waves, $ = Skx + ip): 



F, 



U R = U cos2V, 



0, (Fi) a = F s 



dU R 

dx ' 
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where (• • • ) s denotes the averaging over the microscopic 
spatial oscillations with the period of the order of light 
wavelength. The force F R is rectified gradient force by 
the terminology of [7, 8] as it has the order of magnitude 
of the gradient force, and it is constant-sign on macro- 
scopic spatial scales L ~ 1/Sk much greater than the 
light wavelength A. 

Very similar physical scenarios of gradient force recti- 
fication where some frequency components of the optical 
field induce nonuniform redistribution of atoms over the 
quantum states and other frequency components form 
gradient forces (dependent on atomic state) can be orga- 
nized for polychromatic fields of the 3D spatial configu- 
ration. In particular, this is demonstrated in the present 
paper. In this scenarios, the presence of the spatial mi- 
crooscillations (of quantum state populations and gradi- 
ent force potentials) with close frequencies or with phase 
shift each other is a necessary condition for the gradient 
force rectification. 

So, the simple model described here and more general 
results of works [7-10] show that RRFs possess the follow- 
ing remarkable properties which allow one to use them 
to create 3D dissipative all-optical traps performing cool- 
ing and trapping of the particles. These forces are sign- 
constant at macroscopic spatial scales L greatly exceed- 
ing the light wavelength A (L 3> A), and in strong fields 
they have an order of magnitude of the induced light 
pressure force (gradient force [l|) and are not saturated 
with increasing field intensity (thus, RRFs can signifi- 
cantly exceed the spontaneous light pressure force). The 
other useful property (which manifests both in strong 
and weak fields and which is not paid due attention to) 
is the controllability of the spatial structure of RRFs. In 
particular, choosing properly the configuration of the in- 
terfering light beams and their parameters it is possible to 
create the purely potential RRF (completely suppressing 
its vortex component). Such a RRF forms the system of 



deep potential wells [11| with a depth exceeding the char- 
acteristic value of the light (Stark) shifts of the atomic 
energy levels hA s by a large factor L/A> 1. 

Moreover, in the general case, RRF contains a dissipa- 
tive component — friction force Q. 



The RRFs properties mentioned guarantee the possi- 
bility of complete overcoming OET constraints and re- 
veal the perspective of their use to create new all-optical 
schemes of cooling and trapping of atoms. 

However, certain solutions significantly depend on the 
structure of atomic transitions. The 3D atom localiza- 
tion 85 Rb in an optical superlattice (induced by RRF) 
was demonstrated by experiment in (l2j ]. This paper 
notes an interesting possibility of applying such meth- 
ods of the atom localization in quantum computing op- 
erations. In [l3T ] the authors investigate the possibility of 
using RRFs for 3D all-optical trapping and cooling of the 
atoms with the quantum transition Ft = 1/2 — > F a = 3/2 
(where Fi, and F a are the total angular momenta in the 
ground and excited states, respectively). 

Based on the performed 3D semiclassical Monte-Carlo 
simulations of the atom motion in the biharmonic field 
(3D configuration), it was shown that in the case consid- 
ered the stable atom localization in the 3D superlattice 
can be combined with their sub-Doppler cooling. 

In [ll-Gji the three-dimensional rectification of a radi- 
ation force in weak and strong non-monochromatic fields 
in the case of the atoms with the quantum transition 
Ft, = — > F a = 1 was investigated, certain all-optical 
schemes of the 3D stable confinement of such atoms were 
suggested. The temperature of the trapped atoms in the 
cases considered in EMU can not exceed the value of 
the order of Doppler cooling limit. 

In the present work the problem of all-optical 3D trap- 
ping and cooling of resonant particles is studied using 
RRFs for the cases of the atoms (ions) with the tripod- 
type configuration of the working levels: for the atoms 
(ions) with the quantum transition Ft, = 1 — > F a = 0. 
The case considered is significantly different from that of 
the atoms with the quantum transition Ft, = — > F a = 1 
due to the degeneracy of the ground state. This im- 
portant factor predetermines the possibility of the sub- 
Doppler cooling. On the other hand (as we shall see 
later) the tripod-type atom is convenient for the theoret- 
ical analysis of the 3D problems. 

It was shown that RRFs, induced by non- 
monochromatic optical fields are able to provide the sta- 
ble deep 3D localization of tripod-type atomic particles 
and efficient control of their temperature. 

In the proposed 3D dissipative optical trap the tem- 
perature of the localized atoms can be varied in a very 
wide range of values including both super-Doppler and 
sub-Doppler temperatures. 

The desired effect (trapping and cooling of atoms) can 
be achieved using the 3D scheme for rectifying a gradient 
force based on employing the partially coherent optical 
field involving (besides the coherent components) compo- 
nents with fluctuating phases (the same idea for the case 
of atoms with the quantum transition Ft, = — > F a = 1 
was considered in [15j ) 

A distinctive peculiarity of the investigation conducted 
is the approximate analytical description of the 3D effects 
of the RRF action on the tripod-type atoms in terms of 
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the kinetic quasi-classical theory of the radiation forces 
based on the Wigner density matrix formalism, allow- 
ing one to correctly take into consideration the quantum 
fluctuations of the radiation forces. Note that in most 
works devoted to RRFs (see [Toj]) the ID problems are 
usually analyzed and RRFs are determined using the con- 
ventional optical Bloch equations (the exception from the 
papers cited is the article [Hj]). 

The paper is organized as follows. Section II describes 
the kinetic quasi-classical model of the mechanical action 
of non-monochromatic light on the tripod-type atoms 
and procedure of averaging equations for the Wigner den- 
sity matrix, resulting in their significant simplification. 

Section III describes in detail the 3D configuration of 
the optical fields which provides the desired effect of the 
stable localization and cooling of atoms. 

In Section IV the Fokker-Plank equation for the 
Wigner distribution function of the atoms /(r, v, t) is ob- 
tained, and based on the analysis of its coefficients from 
the governing parameters different characteristic regimes 
of the mechanical action of optical fields on the tripod- 
type atoms are described. 

In Section V an asymptotic velocity distribution of the 
atoms is found and the Smoluchowsky equation for the 
configuration space distribution function n = n(r,t), de- 
scribing the slow (diffusion) stage of the evolution of an 
atom ensemble is obtained. The possibility of the long- 
term spatial localization of the atoms in deep potential 
wells is shown and the dependence of temperature of the 
trapped atoms on the field parameters is analyzed. Par- 
ticularly, the possibility of the stable long-term trapping 
of atoms both with Super-Doppler and Sub-Doppler tem- 
peratures is demonstrated. 



II. MODEL 

Consider an ensemble of atoms in the light field 

E(r,^)e~ l " ot + c.c., 

with the carrier frequency wo, tuned resonant to \Fj, = 
l,M b = 0,±1) -> \F a = 0,M a = 0) closed (cycling) 
atomic transition, where F a , M a is the full angular mo- 
mentum and its projections for the ground a = b and 
excited a — a states. The field is a superposition of the 
coherent quasi-resonant components (with three different 
frequencies), polarized in mutually perpendicular direc- 
tions and partially coherent (fluctuating) resonant field 
E' with the bandwidth ~ T: 

E(r,i)= E jl (r)e J exp[-iA 3 t}+-E'(r,t), (1) 

j=x,y,z 

where denotes the unit basis vectors of the Cartesian 
coordinate system and Aj is detuning from the resonant 
frequency u>o ■ According to the original conception of the 
effect of the radiation force rectification [7|, |8| (compare 



also with [l5|) one can assume the following hierarchy of 
the characteristic frequencies: 

Aj - ) |A J -Ai|»l^il, 
r»|c/il, J^A 7, (2) 

where I and j ^ I denote the indices x, y or z, Vji = 
dEji/h, Uj = d{e.j ■ E'*)/h are the Rabi frequencies, 

d = \\d\\/y/3, \\d\\ is the reduced dipole transition matrix 
element, k — loq/c is the wave number, s is the ther- 
mal velocity of atoms (characteristic width of the veloc- 
ity distribution of atoms), 7 = 7'/3, 7' is the rate of the 
spontaneous decay of the excited state. The relations be- 
tween the frequencies in the right part of inequality ([2]) 
can be arbitrary. Inequalities © _ © imply that the co- 
herent components of the field are "quasi-resonant" , and 
the fluctuating component E' is "resonant". Therefore, 
(as it will be seen later) the coherent components of the 
field, cx Eji(r)ej, form the spatially non-uniform light 
(Stark) shifts of the atomic energy levels, and the fluctu- 
ating field component E'(r, t) provides incoherent excita- 
tion of the atoms (their redistribution over the quantum 
states). In the scheme of the gradient force rectification 
in the bichromatic field considered in 0, Q this effect is 
achieved due to the presence of the coherent field com- 
ponent with a relatively low frequency detuning. 

The state of the atoms interacting with the optical 
field will be described using the Wigner density matrix 
p(r,v,t) [H, 0. In the quasi-classical limit, when the 
photon momentum is much lower than the characteristic 
atomic momentum distribution width, hk -C ms (where 
m is the atomic mass), and in the interaction represen- 
tation, this density matrix satisfies the quantum kinetic 
equation 0, El 

d P ~ -WrA 1 jdVdp] 

i +w= -W + ^(*55)' (4) 

d__ d_ d_ 
~dt~ dt +y> "dr' 

where hV is the dipole atom-field interaction Hamilto- 
nian, 7 is the relaxation operator that includes the recoil 
effect during spontaneous transition [J HI; and the square 
brackets and braces denote the commutator and anticom- 
mutator, respectively. The second term in the right part 
of Eq.Q takes into account the recoil effect in induced 
transitions. Further, it is convenient for our analysis to 
consider p in the Cartesian representation (compare with 
[Til . [l9| ) i.e. in the representation of the basis wave 
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functions (of intra-atomic motion) for the excited \a) and 
ground states \bi): 



|a} = |0,0), |fcs) = |l,0>, 



\bx) 



V2 



\by) = i — 



V2 



In this representation the matrix elements of the dipole 
moment d are directed along the unit vectors of the 
Cartesian coordinate system, 

(bi\d\a) = e t d, i = x,y, z, 

and the Hamiltonian HV in the rotating-wave approxi- 
mation (|Aj|, r -C wo) has the following form: 

hV = -h [{Vji(r)exp(iAjt) 

+U j (r,t))\b j )(a\+H.c] (5) 

For the density matrix elements in the Cartesian repre- 
sentation the following notations will be used: (bi\p\bi) — 



p u (r,v,t), (a\p\a) = p(r,v,t), (bi\p\bj) = q tj (r, v, t) with 
i 7^ j, (bi\p\a) — pj(r,v, t). Thus, pa, p are the Wigner 
distribution functions for the atoms in the states \bi) and 
| a) , respectively, the functions pi determine the projec- 
tions of the induced dipole moment onto the axes of the 
Cartesian coordinate system, characterizes the coher- 
ence between \bi) and \a) states. 

From Eq. © it follows that the coefficients of the equa- 
tions for the density matrix elements (see Eq.©) os- 
cillate with the frequencies Aj(j = x, y, z) and, more- 
over, contain fluctuating (stochastic) components. Tak- 
ing into account inequalities ©,© we average succes- 
sively Eq.© — at first, over the high-frequency oscilla- 
tions with the frequencies Aj (as is usually done in the 
RRFs theory (see @, 0, UJ), then, over the fluctuations 
of the field E'. At the first stage one obtains the system 
of equations (using the same notation for the averaged 
quantities): 



dt +l1 



iAi(r) 



E QH u j 



Ai 



(6) 



— - iAij(r) ) qtj + ijS, 



. , E qit = + (Pi U j ~ U iPj) + % E 

l=x,y,z l=x,y,z 



{pxUf - c.c.) + B 



i j > 



(7) 



dt 



m 



E 

=x,y,z 



+ C.C. 



E 

=x,y,z 



dqa 
dv 



D, 



f E 



(8) 



where /(r, v, t) — Sp(p) = p + J2i Pa i s tne Wigner par- 
ticle distribution function (DF) in the phase space (r, v), 
9« = Pa — p denotes the densities of the distribution of 
the population difference, 



Ai(r) = - 



2|Vh(r) 
A, 



l^iWI 5 



fhk 


n 




V m 







are the effective spatially nonuniform (in case of the non- 
homogeneous fields) detunings due to the light-induced 
Stark shifts, A I3 (r) = A,(r) - Aj(r), 



(9) 



is the operator describing the recoil effect at spontaneous 
transitions t 2l|. The terms Ai and Bij in the right part 
of Eqs. ©, © describe the influence of the recoil effect 
on the induced transitions and coherence between \bi) 
and | a) states. Explicit expressions for them are given 
in Appendix, as well as some comments on the effects 
described in Eqs.©-©. 



I 

Next, Eqs. ©, © are averaged over the fluctuations of 
the field E' [22| . These equations are the system of mul- 
tiplicative stochastic linear equations and the standard 
procedure described in [23| was used to average them. 
In the problem under consideration it is based on the 
expansion of the solution of Eqs. ©, ©, © in terms 
of £ <C 1, which is proportional to the correlation time 
Tc-T- 1 : 

\Uj\r c , |Aj-|r c , ksr c , jt c < ( < 1. 

The smallness of the parameter ( is due to the original as- 
sumptions: the right part of inequality ©, i.e. due to the 
conditions which determine the value of the bandwidth 
r of the field E'. Moreover, one makes the following 
additional assumptions on the properties of the random 
processes Uj, j — x, y, z: 

((Uj)) = 0, «tf J (r ) t)E7i(r',i + T)>) = C i (|r| ) r ) r')* Ji) (10) 

where the double angular brackets denote averaging over 
the fluctuations. 

Thus, suppose that Uj is treated as a stationary ran- 
dom process with zero-mean, and the E' components 
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with different poiarization fluctuate independentiy. As 
a resuit, the following closed system of the reduced equa- 
tions describing the kinetics of the tripod-type atoms in 
the light field was obtained (without changing notations 
for the averaged DF): 



J t +^(*)+l 



df 9 



= jf + B h i = x,y,z, (11) 




(12) 



The system of four equations (|TT|) — (|T2|) is much more 
simple than the system of sixteen ones ([S])-© for the ele- 
ments of the density matrix. Such a simplification results 
from the original assumptions on the relations between 
the parameters of the optical field and the properties of 
the field component correlators: i.e. from inequalities 
(J2) (O and Eqs. (fTU)) . Following the generally accepted 
approach of the quasi-classical theory of the mechani- 
cal action of light on atoms [l|, , later (in the Section 
IV) one obtains from Eqs. (fTTj) . (IT2|) the Fokker-Planck 
equation (FPE) for DF /(r, v, t); however, to obtain the 
explicit FPE representation it is necessary to specify the 
spatial configuration of the optical fields. In the next 
Section, in accordance with the aims of this work, such 
a field configuration will be taken where the 3D kinetics 
of the atoms is completely determined by RRFs. 



s{Q} = £ 



i,3,j' 



A i 



(13) 



where q t = fa)), Q i — ((Qu)), Qu is determined by 
the expression (|A.1[) . Vj = v • e^, and Ri, Fj, Ajj, are 
determined by the correlators 



III. OPTICAL FIELD CONFIGURATION 

Let the optical field be formed by a special superposi- 
tion of the plane light waves, where the Rabi frequencies 
Vn and Ui (i — x, y, z) are the following: 

Va(r) = y [exp(iq, • r + %) + exp(iq^r)] , (18) 



Jfc(r) = 2Ke / ((Ui(r,t)U?(r,t + r)))dr, (14) 



Ui(r) = £ Ui a (r) exp (iip a (t)) , 



a=l 



(19) 



Fi = ~2hSm / ((VUi(v,t)U*{r,t + r)))dr, (15) 



A 1 



2m 2 



5Re 



dU*(r,t) dUi(r,t + T) 



dr-j 



dry 



dr. 



(16) 



The coefficients i?.;(r) have the meaning of the rates of 
the transitions between the low- lying \bi) and excited \a) 
atomic states induced by the field E', the value Fj has the 
dimension of force and is proportional to the energy flux 
density of the E' field component, polarized along e*, the 
coefficients A*-, v have the dimension of velocity diffusion 
coefficient and the order of magnitude (hk/m) 2 Ri. The 
last term in the right part of equation (fTTj) Bi = —i((Ba)) 
(where Bu is determined by Eq. (|A.4[) and has the fol- 
lowing explicit representation 



B, = 



_L V — 

4to i9v 

| 1 dQ t 
2m <9v 



(Qj - 2qj)V 



(Qj 



,h\V a 



(17) 



Un(v) = — [exp(ikii • r) + aexp(ik' ll ■ r)] 



U i3 (r) 
U t2 (r) 



>a'U 



[exp(«k i3 • r + <i) + exp(ik- 3 • r) 



t/A(r), Uu(v) = U; 3 (r), 



(20) 



where the index i — x, y, z, rji, and £i are the fixed phase 
shifts, the positive parameters a, a' <C 1, U and Vi are the 
real amplitudes (which in the considered model do not de- 
pend on r), <pia(t) are independently fluctuating phases 
(with delta-correlated zero- mean derivatives), which de- 
termine the correlators of the E' components by the re- 
lations 



((expi((pj a (t) - (fip(t - 



fy<J a/3 exp(-r|T|). (21) 



Thus, the field E' is described by the so-called phase- 
diffusion model [It], [U, HH and has a Lorentzian spectral 
profile J(oS) with bandwidth V: 

r 

J(u>) oc 



( w -^ ) 2 + r 2 - 



Here, all the assumptions on the properties of the random 
processes Uj are implemented (see Eqs. (ITU)) ). 

It is worth noting that the representation (fT9"|) for 
Ui(r, t) is valid only if the coherence length l c = ct c = c/T 
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FIG. 1. The wave vectors of the light field components po- 
larized along et and determining the Rabi frequencies in Eqs. 
(|18H20p . the ordered index triplets (ijl) are: (xyz), (zxy) and 
(yzx), Si and 8[ being small angular detunings. The bold 
arrows are the wave vectors corresponding to the coherent 
field components. The vectors ki2 = — kii, k^ 2 = — k^ij 
= — k;3, ki 4 = — ki 3 determining Ui2(v) and Ui4.(v) are 
not presented in the Figure. Thus, each component E[ei of 
the field E' is a sum of four pairs of the counter-propagating 
waves. 



is much greater than the characteristic atomic cloud size 
D (see detailed discussion in [24j |) 



l c > D. 



(22) 



Specifying the configuration of the optical field, one 
determines the wave vectors in Eqs. (fT8l) as follows (see 
Fig.l): 

kii = k(ej cos (3 + ei sin /?), 

Ki = fc( _e j C0S P + e i sm /3); 
q 4 = q l (e J cos Pi + ei sin ft), 
q- = qi(—ej cos /3, + e; sin/?,), 

k i3 = fc(ej cos/3^ + ej sin /?,'), 

l4=fc(-ejcos i 8 i ' + ej8in/3j), (23) 

where /? is an angle determining the dominant direc- 
tions of the propagation of the waves polarized along e,, 
Pi = fi+5i, P[ = P+S[, Si and S[ are small angular detun- 
ings (\Si\, | <5- 1 <C 1), and the ordered index combinations 
(ijl) are: (xyz), (zxy), (yzx). The angular detunings 5, 
and are free parameters allowing one to control the 
macroscopic spatial structure of the field (compare with 
0), we choose them in such a way that the following 
conditions are satisfied for the difference of the wave vec- 
tors: 



Aqi = q, - q ■ = k i: 



k' — k' 



2n 



V4, 



(24) 



AKj = Aqj-Ak 4 = 2fc«ift cos f)ej = — sgn(ftaj) cos/3ej-, 

Li 

where a t = (Aj/wo) < 1, Ak 4 = k,i - k^ = k- 2 - k i2 , 
the vectors k l2 = -k,i, k- 2 = -k- l5 k i4 = -k i3 , 



k- 4 = — k- 3 determine Ui2(r) and J7a( r ) in the super- 
position (fTT))) (in accordance with Eqs. (|2T)10 . the pa- 
rameters ft = (1 — 5f/2cti — Si taxi p) determine the 
macroscopic spatial scales of the problem Li (i.e. as 
one will see, the period of the spatial RRFs modulation) : 
Li = 7rc/|A,*ft| 3> A = 2ir/k. Further, we will restrict 
ourselves to the case of the configurations of the optical 
fields and values of the angles /3 close to tt/2 or 0, for 
which the following relations hold 



Li 



L, i = x,y,z, Pi = |sin2/J| < 1. (25) 



Then, given the values Aj and L » ) and not very 
small values of P (tan/3 ^S> \/|a!ift|), the necessary (to 
fulfill Eq.(f2"3]) and the first one from Eqs. (US])) values of 
the governing parameters Si and 5[ are described by the 
simple formulae: 



(1 - 
tan p 



a»ft 
tan /3 ' 



(26) 



where |ft| = 7rc/|A;|L, and the sign ft may be taken 
arbitrarily. 

So the following conclusions can be drawn: each group 
of light waves polarized along i-axis can produce an in- 
terference pattern along j-axis (where index j = y for 
i = x; j — x for i — z; j — z for i = y). Spa- 
tial frequencies of this interference patterns Aq ; and 
Aki are close to each other (due to smallness of angu- 
lar detunings Si, 5^) and determine microscopic spatial 
scale of the problem Am — A/2 cos p ~ l/|Akj|, l/|Aq,| 
while their difference AK^ determines macroscopic spa- 
tial scale: L ~ (l/|AKj|) > \ M . 

In the field configuration considered the explicit ex- 
pressions for the transition rates Ri(r) and the term S(Q) 
in the right part of equation (fT2")) for DF f(r, v, t) are the 
following 

U 2 

Ri(r) = R(l + aiPi(r)), R=—(l + a 2 + 2a'), 



S{Q} 



[cos( Ak, • r) + b cos( Aq, • r + ft )] , (27) 




d 2 d 
cos 2 P— -j +sin 2 /3— T 
OVj ovf 



(28) 



where Oj = 2(a+a')/(l + a 2 +a'), b = a' /a and in Eq.(| 
the ordered index triplets (ijl) are: (xyz), (zxy), (yzx). 
The most important peculiarity of the field configuration 
determined by Eqs. (IT8)) -(l20 | is vanishing of the forces of 
the resonant light pressure, conditioned only by the fluc- 
tuating field E'. This results from the accurate mutual 
compensation of independent (due to Eqs. (|2"Tj) ) contri- 
butions into Fj (see Eqs.(fT5]), proportional to the energy 
flux density of the field components with the Rabi fre- 
quencies Ui a , a = 1 — 4: 

F< = 0. (29) 
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Eqs. (1291) hold in case when Eqs. ([20]) are satisfied. 

Now one can clearly see the final distribution of the 
"roles" of the coherent Ei and partially coherent E' fields 
in the considered model of the mechanical action of light 
on the tripod- type atoms, see Eqs. (jTTJ) , (fT2")) . The fluc- 
tuating field E' is responsible for incoherent mixing of 
the atomic states (or in other words for the redistri- 
bution of the atoms over the quantum states), and the 
quasi-resonant coherent field induces the effective poten- 
tials (due to the light-induced Stark shift of the energy 
levels), which determine the motion of the atoms: the 
unexcited atoms (in the state \bi)) move in the field of 
the gradient force with the potential (oscillating with a 
period of the order of the light wave length A) 



n\v t ir\' 

A,- 



x,y,z, 



the excited atoms (in the state |o)) in the field of the 
gradient force with the potential 



y, h\Vg(r)\< 

=x,y,z 



A, 



In order to better understand the above-described 
physical picture, it is sufficient to pay attention on the 
fact that Eqs. (11), (12) (at condition Eq. (29)) are equiv- 
alent to the following equations for Wigner distribu- 
tion functions pa and p for the atoms in the states \bi) 
(i = x, y, z) and \a) {p=(f -J2 <fe)/4, Pu = P + qu) 



dpu 
dt 



dp 
~dl 



VhlVnl 2 d Pi .. 
A,- 9v 



jp-R i (r)(p ii -p)+S i , i = x,y,z, 



E 



VftlV^il 2 \ dp 
A t dv 



where Si and S are second-order terms in the quasi- 
classicality parameter hk/ms responsible for the veloc- 
ity diffusion (conditional on spontaneous transitions and 
fluctuating field E' influence). These equations clearly 
demonstrate four-potential nature of atomic kinetics and 
a factor of the incoherent redistribution of population of 
quantum states (relevant to these four potentials). The 
mean gradient force is obviously, a sum of gradient forces 
weighted by the probabilities of occupation IT = Pi/ f 
(i = x,y, z), II = pi f of atom states. 

Therefore, the physical model proposed is a general- 
ization of the simple model of the mechanical action of 
bichromatic standing wave on two-level atom, which is 
described in the Introduction. The only difference is the 
number of atomic states, more complicated spatial struc- 
ture of gradient force potentials and much more complex 
nature of atomic population redistribution among quan- 
tum states. 

Taking into consideration this analogy it is possible to 
assume that the necessary condition of the gradient force 



rectification in our problem (cf [15]) is the spatial modu- 
lation (with the period ~ A) of the transition rates Ri(r), 
and hence the spatial modulation of the relative popula- 
tion differences fa = qt/f, which is possible only if the 
field E' has mutually interfering components, i.e. when 
a, a' ^ 0. Though we consider the weakly nonuniform 
fluctuating field E', a, a' <C 1, the effect of the RRF ac- 
tion on the atoms can be strong, particularly, due to the 
ability of RRF to retain the sign on macroscopic spatial 
scales L 2> A. In the next section we will see in more 
detail how the rectified gradient force is calculated and 
how the RRFs act on the tripod- type atoms. 



IV. FOKKER-PLANCK EQUATION (FPE) 

One obtains first of all the kinetic equation FPE, de- 
scribing the DF / evolution at times 



Ri 



R 

cos j3 



(30) 



Use is made of the procedure of the adiabatic elimina- 
tion of the internal degrees of freedom of the atom which 
is well-developed in the theory of the resonant light pres- 
sure [IHl] , i.e. in our case, the elimination of the variables 
qi, included in Eqs. (|TTj) , (fT2"|) . In fact, it implies the sep- 
aration of the fast processes of the redistribution of the 
tripod-type atom over the quantum states (occurring at 
times t ~ i? _1 ) from the slow processes associated with 
the translational motion (occurring at times t ^> i? _1 ). 

One should take into account the smallness of the pa- 
rameter ai — 2(a + a')/(l + a 2 + a') (according to the 
assumption on the weak non- homogeneity of the field E', 
at a, a' <C 1, see Eqs. (|2"0"|) . (|2"7|)) and quasi-classicality 
parameter e = (hk/ms) <C 1. Moreover, we restrict our- 
selves to the case of the slow atoms: 



2 9 



(31) 



where v c ~ R\/k has the meaning of a so-called cap- 
ture velocity, i.e. of the characteristic velocity, which 
determines the area of the most efficient RRFs action 
in the problem under consideration: |v| <C v c . It can 
be shown that outside this region the rectified gradient 
force (RGF) and coefficient of the light-induced friction 
fall rapidly proportional to (y c /v) 2 . The physical mean- 
ing of formula v c ~ R\/k = R/kcos/3 is very simple: v c 
is a typical velocity at which the atom travels a distance 
of the order of microscopic spatial scale A m = A/2 cos (3 



during the characteristic time ~ R , i.e. 



XmR] 



the slow atoms travels a small distance (compared with 
the Am) during this time. That is why v c value can be 
adjusted by geometric parameter changing: cos/3 (this 
effect is described in the paper [20], Section 2.3 for two- 
level atoms and bichromatic field). It is worth noticing 
that in the case of the rectification schemes considered in 
[S EE v c is determined only by the spontaneous re- 
laxation rate and does not depend on the field intensity: 
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v c ~ l/k- The quasi-stationary solution of Eqs. (|TT|) . de- 
termining the relation between and / at t 3> 7 _1 , -R -1 
is presented as an expansion over the small parameters 
of the problem (e, a%, (v/v c )): 



Qi 



-(0) t i -(1) r 

8i '/ + ?< / 



(32) 



where q^ are linear operators. 

Using linear approximation over the velocity v and pa- 
rameters e, a\ (i.e. taking into account in expansion 
(l32l) onlv the dominant terms, proportional to the small 
parameters to the power not exceeding the first power) 
the following expression for qi was obtained: 



/ 



4 X + 3 



E 

a=0 



f, 



(33) 



Qio = 



aif 
4 X + 3 



[2(2 X + l)P t ( T ) - P t (r) - Pfr)] 



Qn 



m = n 



aif 



(4 X + 3) 2 i? 



[2(6 X 2 +8 X + 3)P,(r) 



^(4 X 2 +8x + 3)(F i (r)+F,(r))], 
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mdv (4 X + 3) 2 i? 



?i3 



N 



(4 X 2 + 5 X + 2 
-(2 X + 1)V 



ft|Vn(r)| 2 , 



1 d 

7 (4 X + 3)S' 



fi|^-i(r)P 
A,- 



where I =^ j =/= i, P^(r) = v ■ VP;(r), X = P/7. The 
second term in Eq. (|33f describes the effect of the non- 
homogeneous distribution of the atoms over the quan- 
tum states due to the spatial modulation (with the pe- 
riod ~ A) of the transition rates Ri(r) and light shifts 
|Vii(r)|/Aj. It is important that non-local effects (the 
effects on non-adiabaticity in the atom response to the 
action of external fields) are taken into account, since the 
terms qi a f (a = 1 — 3) in Eq. (l33l) are determined by 
the field intensity gradients and are connected with the 
atom motion and recoil effect. 

Substituting expression (|3"31 for qi in Eq. (IT2"j) , re- 
taining terms up to the second order in the parameter 
e = hk/ms and averaging the equation for DF over the 
small-scale spatial oscillations with a period of the or- 
der of the light wave-length ~ A, one obtains after some 
transformations the following FPE for the Wigner distri- 
bution function averaged over the microoscillations, /: 



df 
dt 



d 



md~v 



(F R + F m )f = D 



dvf 
(34) 



where F# does not depend on the atom velocity and has 
the meaning of the rectified gradient force (RGF), 



E 

i=x,y,z 



RV|V5i(r) 



-&o(r) 



F#i is linear-in- velocity RRF (the retarded gradient force 
by the terminology of [1] , which appearance is associated 
with hysteresis in the response of the moving atom on 
external field): 



'Rl 



E 



fiV|Vii(r) 
A, 



-&i(r,v) 



D, Dm are the diffusion coefficients (determined by the 
terms q^f and qaf in the expansion (33) and also by 
terms of the right side of the Eq. (12)). Note, that in our 
approximation (|33|) we neglect small additives to RRFs, 
which have the order of magnitude a\ and {v/v c ) 2 ai. 
The coefficients of Eq. ([34|) do not contain components 
oscillating (in space) on microscopic spatial scales and 
can change only on macroscopic scales > L. This is ex- 
plained by the fact that relative differences of population 
Qi = lit f (See Eqs. (27), (33)) contain oscillating compo- 
nents with spatial frequencies close to spatial frequencies 
of gradient force oscillations ex V|Vii(r)| 2 . The condition 
of the correctness of the implemented procedure of av- 
eraging over the small-scale spatial oscillation (which is 
quite similar to the case of the two- level atoms [2(|) is 
the following: 



ms 



ftl^WI 2 
A,(4 X + 3)' 



(35) 



and implies that the effective temperature T (in energy 
units) of the atoms is considerably higher than the depth 
of the macroscopic potential wells produced by the os- 
cillating (with a period of ~ A) gradient force with the 
potential 



U g (r) = 



V h\Vu(r)\ 2 



In the situations considered here the condition is al- 
ways satisfied and the rapidly oscillating component / of 
DF is a small correction to /: (/// ~ ai <C 1). 

Given all the detunings A, > 0, the parameters < 
(i = x, y, z) and the Stark shifts of the energy levels in- 
duced by the coherent fields with mutually orthogonal 
polarization directions arc equal, i.c the the following re- 
lations hold: 



A x 



I 



v 2 



(36) 



where / = I x is the intensity of the light waves forming 
the field component E x ie x , I s = hujok 2r y /6ir is the in- 
tensity of the optical radiation saturating the quantum 
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transition. In this case the explicit expressions for the 
forces included in Eq. (|34l) . have the following compact 
form: 

F R = -VU R + F RQ , Fju = -m *( r i)viei, (37) 



U R = U ^2 cos (j^ r i - Vij 



u = hj 



kLi aiG 
~2n (1 + b) 



Mx), 




' RO 



(38) 



Ci)A(x) J2 *> ( 39 ) 



i—x,y,z 



2(6x 2 + 8x + 3) aiGcos^ 
(4x + 3) 3 X (l + b) ' 



2tt 

6cos(?7i - Ci) + cos ( —r.i - 7/1 



. (40) 



where r l = (r • e,), u,; = (v • e, ; ), L x = Lj cos/?, A(x) = 
(2 X + l/(4 X + 3) 2 . 

One can see from Eqs. ([3"? ]) -([3"9" ]l that RGF F R con- 
tains in the general case (771 7^ Ci) the component Fro 
which does not depend on the spatial coordinates r. 
This is the manifestation of the total rectification effect 
appearing due to a special choice of the wave vectors 
(see Eq. (|2"4")l and Fig. 1). Another component of RGF, 
—\7Ur, gives rise to an optical superlattice, a system of 
periodically (with a period L\ S> A) distributed poten- 
tial wells for the atoms. Though the force F^o is likely to 
have special practical applications; here, one is interested 
in a certain problem of the 3D localization and cooling 
of the atoms, and, therefore, should restrict oneself to a 
special case of choosing the phases of the fields r\\ and £1 
and the parameter b — a' /a in Eq. (|4T)j) : 



m = Ci 



b> 1. 



(41) 



Thus, Frq = 0, and F^i is a linear-in-velocity friction 
force, determined by the spatially nonuniform strictly 
positive friction coefficients 



Xi = x(r») > 0. 



(42) 



Interestingly, in this case expressions (|3"T)) - (|40I) coin- 
cide with the expressions for RRFs, obtained for the 
tripod-type ion by another method in (26[. 

Thus, under the conditions considered RRFs induce 
the 3D dissipative cubic superlattice, capable of simul- 
taneously cooling and trapping atoms. However, the ki- 
netics of the atom cooling and trapping considerably de- 
pends on the values of the velocity diffusion coefficients 
D and Drj (determining the process of the atom heating 



competing with the atom cooling) 
explicit expressions were obtained: 



Here, the following 



D = D 



si 



. hk 

D M , ~ — 



IX 



m J (4x + 3)2' 



D 1 c - 



hk 



mj 7 (4x + 3)4 ' 



x(4x + 2) 



D 



16x 2 + Hx + 3 



hk\ 2 8x^ 

ray 7 (4x + 3) 3 x 



G 



D, 



(43) 



where j = x,y,z, D s \ is the diffusion coefficient condi- 
tioned by the recoil at spontaneous transitions, D\ is the 
coefficient of the induced diffusion connected with the 
action (on the atoms) of the fluctuating field E'[27j ( its 
order of magnitude coincides with that of the diffusion 
coefficient in the case of the two-level atoms in the field of 
the standing wave with the fluctuating phase [13), D R is 
the coefficient of the induced diffusion connected with the 
quantum fluctuations of the gradient force. When writ- 
ing Eqs. (43) we neglect small corrections of the diffusion 
coefficient ~ a 2 , (s/v c )a\. 

The distinct feature of the considered 3D scheme of the 
mechanical action of the light on the tripod-type parti- 
cles is the simultaneous presence of several independent 
governing parameters: x, G, L±, cos/3, 01, b, Ci, Vi- The 
specific choice of these parameters determines the values 
of the force and diffusion coefficients of FPE, and, con- 
sequently, various regimes of the atom localization and 
cooling. Now we shall describe some characteristic situ- 
ations, considering Eq. (|4*Tj) to be true and using Eqs. 

(133, (US), m, and m. 



A. Strong fields: G > 1, \ > 1 

In this case RGF Fr ~ F s Gai/x and can considerably 
exceed the maximum value of the so-called spontaneous 
light pressure force F s — Hkj/2 [1], if G ^> x/ a i- The 
capture velocity v c can exceed the characteristic value 
v c0 = j/k by a large factor (x/ cos/3) 1. The friction 
coefficient can exceed the characteristic value ljr: h ~ 
ujrcii cos/3G/x 2 ^> (*)r, if G 3> x 2 /°i cos/3. The depth of 
the macroscopic potential wells AU = 2Uo considerably 
exceeds the depth of the microscopic potential wells U g : 
AU ~ U g (Ga 1 L/\x) > U g , in case (flti-L/A) > x/G. 

Thus, high values of RGF and capture velocity v c can 
be achieved as well as super- fast cooling rates, if G ^> x 2 - 
However, in this regime the atom heating rate is also 
high. Really, comparing the velocity diffusion coefficient 
D = D + Dr (determining the heating rate) with the 
characteristic value of the diffusion coefficient in a trav- 
eling plane wave [IHi, D = (hk/m) 2 7, one obtains: 



D = Do 



91 

»X 



X » D . 



10 



B. Strong coherent field and not very strong 
partially coherent field: G > 1, x ^ 1 

One has the following relations: Fr ~ F s Ga\ 3> F s (if 
Goi > 1), AU ~ U g (L ai /X) > £7 S (if L > A/ai), x - 
w R Gai cos/3/x > (if G > x/aiCOS/3), w c = w c0 (if 
X ~ 1 cos/3 w 1) and v c ^> w c o (if X ~ 1 and cos /3 <C 1), 

D ~ Z> G 2 /x » A>. 

This implies that in this case it is also possible to reach 
super-high values of RGF, high values of the capture ve- 
locity v c (but only by means of the proper choice of the 
geometric factor: cos/3 <C 1), super-high cooling rates 
and intensive diffusion in the velocity space. The rela- 
tively small intensities of the fluctuating field E' make it 
significantly different from the case A. 

A similar regime of the mechanical action of the co- 
herent bichromatic light on the two-level atoms is de- 
scribed, for example, in: [M, l20j|. however, these descrip- 
tions cannot be applied directly for the case of the 3D 
configuration (considered here) of the optical fields and 
tripod-type atoms. 



C. Weak fields: G ~ \ < 1 

There exist the following estimates: Fr ~ F s a\G, 
AU ~ U g a\L/\ v c — v C oxl cos /3, x ~ cj^ai cos/3, 
D ~ D G. 

The most significant peculiarity of this case is that the 
friction coefficient does not change with the decreasing 
intensity of the fields, but the coefficient of the velocity 
diffusion decreases. In other words, it appears possible 
to considerably suppress the radiation force fluctuations 
when maintaining the high cooling rate. A similar phys- 
ical situation appears in the famous ID model of the po- 
larization gradient cooling (see [28| and references) and 
it is a prerequisite of achieving the Sub-Doppler temper- 
atures. 



V. COOLING, SPATIAL DIFFUSION AND 
LOCALIZATION OF ATOMS. 

An important feature of the Brownian motion of 
atoms, described by FPE (|34|) for x > is the over- 
damped character of this motion: 



— ~ A r Li, 
x 



(44) 



where A r has the meaning of the effective free-path length 
of the atoms in viscous "fluids" of photons. A similar si- 
tuation occurs in OM (5j], if L\ denotes its dimensions. In 
the problem considered x > ujRa\ (see section IV), there- 
fore, assumin g k ~ 10 5 cm" 1 , 7 ~ 10~ 8 s^ 1 , m ~ 100 
amu, s ~ \fhrjjrn, a\ = 0.2, one has the following: 
A r < 10~ 3 cm. Thus, condition (|44l) of the strong dissi- 
pativity of the optical superlattice is well satisfied, if the 
period is L\ > 0.01 cm. 

Assuming inequality (|4"4"f to be satisfied let us con- 
sider the evolution of the atom ensemble at times t > 
Td ~ x~ 1 (Li/\ r ) 2 ^> T ~ x^ 1 , where (according to 



the general properties of FPE [20|) To is the time nec- 
essary for the local quasi-stationary velocity distribution 
to be achieved, Td is the characteristic time of developing 
a much slower (diffusion) process of changing the atom 
density, n(r,t), in the configuration space (r). 

Then, it appears possible to present the solution of 
FPE (|34l) in the following form (using instead of the vari- 
ables (r, v) the variables (r,c = v— u(r,t)), where u(r, t) 
is the macroscopic (directed) velocity of the atom mo- 
tion) : 



/ = n/ + h 



(h)c (Ac)=0, A«n/ 



(45) 



(46) 



where the angular brackets (• • • ) c denote the integration 
over the velocities of the chaotic motion c, 



! m 



exp 



2T, : 



(47) 



Tj = T(ri) = m{D + Dr) / ' x{ti) — T Ti has the meaning 
of the effective local atom temperature (in energy units), 
characterizing the chaotic atom motion along the axis 
i (i = x,y,z), To is the space averaged effective atom 
temperature, % — T(rj) is the function describing the 
effective temperature dependence Tj on the correspond- 
ing spatial coordinate Tj. Explicit expressions for them 
are obtained from Eqs. (|4"0")l . (|4U)) . taking into account 
Eq. flUJ): 



T = 



m(D + Dr) 
x 



( 2-KTj \ 
L X ) 



h = Vb 2 - 1, 



X Q 



2(6x 2 + 8x + 3)ai6iGcos ; 
(4x + 3) 3 x(l + &) 



-UJR. 



(48) 

In the considered approximation, A r <C L\ the density 
n(r,t) and macroscopic velocity u(r, t) are governed by 
the equations (continuity equation and the force balance 
equation) : 



— + drv(mi) = 0, 
at 



d(nT z ) dU h 



dn 



di 



mXiUi 



0, 



(49) 



(50) 



where x; t — x(ri), Ui = (e^ • u). It can easily be proved 
that /o satisfies the equation (c$ = (e, • c)): 



d 



T\_d_ 

m dci 



and a small correction fx to n/o in the first approxima- 
tion satisfies the equation 
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9 x £ |v m-CiCj duj ^ / mcj \ f dui dT, t \ ^ 5Tj / mcj \ \ 

Lh = /on + ? {— - l ) [W< + ^moVj + ? 2^ I— - 3 J J ■ (51) 



Expanding the solution of Eq. (|51|) into the eigenfunc- 
tions of the Fokker-Plank operator L and taking into ac- 
count Eq. (|49|) . (|50p one can obtain a quite strict estimate 



A 



dr^o Tidri>c 



< 1. 



Thus, it is seen from Eq. gZ]), (08]), that at < > xT 1 
and condition (144)) the atom distribution over the ve- 
locities (of the chaotic motion) is, generally speaking, 
spatially nonuniform and anisotropic (with the excep- 
tion of the points r' located on the bisectors of the unit 
cells of the optical superlattice, in which cos(27rr^/Ti) — 
cos(27rr^/Li)). However, within the limit ft » 1, T « 1 
and T ~ To, DF /o is the equilibrium Maxwell distribu- 
tion with the temperature To. 

Further, combining Eqs. (j49| ~(|50 ]l one obtains the 
Smoluchowski equation (SE), describing the atom diffu- 
sion in the field RGF 



dn 
~dt 



^ dr. 

i—x,y,z 



d 1 



T 



dn 
dn 



(52) 



where Ui = —Uocos(2nri/Li). 

The difference of SE (JS2]) from the classical SE (well- 
known in the theory of the Brownian motion [30( ) is the 
spatial non-uniformity of the friction coefficients xi = 
x(rj) and effective temperatures T = T(n). 

Taking into account the form of the obtained SE (|52|) 
a very significant and interesting circumstance appears: 
in the model considered the macroscopic (diffusion) par- 
ticle motions along the axes of the Cartesian coordinate 
system can entirely be separated (i.e. they can be inde- 
pendent from each other). This is revealed in the station- 
ary solution of SE (|5"2"j) being the product of the quasi- 
Boltzmann distributions: 



n 

i=x,y,z 



T(n) 



exp 



To 



(53) 



where ip^ - are the constants, and the effective potential 
U e ff is determined by the expression 



Ueff(r) 



Uo 
h 



^6 cos 


27IT 


1 2 


2nr 




[ Li _ 


--cos 2 




) 






Ui . 





(54) 



and describes the potential wells with the depth AU = 
2U b/h. 

The condition of the deep particle localization has the 
form: 



r = exp 



AU 



To 



> 1. 



(55) 



A prerequisite of satisfying condition (f55"j) is the pres- 
ence of the large multiplier Ti/A > 10 4 (at L\ ~ 1 cm) 
in the expression for Uq (f3"5| . However, stationary lo- 
calized solutions of SE ([52)) are strong idealizations due 
to the irreversible diffusion particle escape from the in- 
terception area of the real laser beams with the finite 
transverse dimensions. Consider a more realistic model 
allowing one to take this effect into account and to show 
the existence of the spatially localized quasi-stationary 
(long-lived) states of the atom ensemble. Considering, 
for cirtainty, a single cubic cell of the superlattice (de- 
scribed by equation |38|) at Wi = (i = it) with the centre 
in the point r = as the dissipative optical trap (DOT) 
for the atoms, we assume that the absorbing boundary 
condition is satisfied: = [2!| on the boundaries £ 
of this cell. It means that the particles are removed from 
the DOT as soon as they reach its boundaries. 

Then, the long-lived localized states of the atom en- 
semble (corresponding to the lowest (the slowest) dif- 
fusion mode), can still be presented as (1531) . assuming 
that the multipliers ^ = ^(r,) exp[— t/n], where ^(r) 
is the eigenfunction, corresponding to the lowest eigen- 
value Ai = 1/ti of the Sturm-Liouville boundary problem 
(SLP) 



H 



1 $ 



#(r) =0, * ( I = 



H 



dr V rtiK dr 



exp 



Ueff(r) 



To 



(56) 



Using (with slight modifications) the results of the 
analysis of a similar SLP (in the case of T ^> 1), made in 
[26| , one obtains the following expression for the lifetime 
r = ri/3 of the atoms in DOT 



Td 



3ttT(0) lnT 



exp 



AU 



To 



(57) 



where Td — (b+ l)Llmxo/Tobi is the characteristic dif- 
fusion time in DOT in the absence of RGF. Thus, condi- 
tion (I55j) is simultaneously the condition of the long-term 
atom confinement in DOT (for a time considerably longer 
than the diffusion time r<j). Note here that at lnT 3> 1, 
^(r) sa const almost everywhere except for the narrow 
regions (with the width of ~ Ti/VlnT) near the bound- 
aries £ of the trap. 

Consider a number of characteristic properties of the 
spatially localized states of the atom ensemble resulting 
from the constructed model. 
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A limitation for the effective mean temperature To of 
the localized particles results from Eqs. (j4*3"|). 



T > 



hj b+1 
2a ibi cos 



G. 



(58) 



This inequality, on the one hand, pre-determines sat- 
isfying important condition (|35j) due to the smallness of 
the parameter a\\ (2ai/3) <C 1. On the other hand, 
from inequality (|58l) it follows that reaching Sub-Doppler 
temperatures Tq <C Td ~ ^7 is possible only at rather 
small values of the coherent field intensities (i.e. at 
G < {2a\bi/b + 1) -C 1). It seems quite natural since 
in the limit Gel and X •C 1: 



T = BT D \G 



3r; 

G 



> 2V3B X T D = T„ 



(59) 



where Td = tvy is the characteristic temperature, de- 
termining the Doppler limit of the atom cooling [lJ-Q , 
B = (b + l)/2ai&i cos/3 and the equality is achieved at 
G = s/3 X - Thus, the Sub-Doppler cooling regime can 
be implemented at X < a\b\ cos/3/(& + l)-\/3- On the 
other hand, the parameter \ (at G = \/3x) should be 
limited by lower bound: X > w i?.V / 3^i/7ai(6 +1), to sa- 
tisfy condition (|3"Tj) . s 2 < v 2 = R 2 /k 2 cos 2 /?. Regarding 
this circumstance and Eq. (|59p one obtains the estimate 
of the lower limit T' m of the effective mean temperature 
of the described localized states of the atom ensemble 



T 

m 

Td" 



3ujr 
7a 2 



(60) 



It follows from Eq. ^ that TJT D < 1, if {u> R /i) < 



Note, that the typical values of ujr/j are < 10 



-3 



Reaching Sub-Doppler temperatures easily corresponds 
to the condition of the long-term deep spatial atom lo- 
calization (|55|) making the proper choice of the period L\ 
of the optical super-lattice (induced by RGF): 



lnT 



4a 2 6cos/3 Li 
9(6+1) 



(61) 



Now let G 3> 1 and X 0.1. This case corresponds to 
Super-Doppler temperatures of the localized atoms: 



T = T D BG > T, 



D ■ 



(62) 



Moreover, the condition of the long-term deep localiza- 
tion is also expressed by inequality (l6TT) . It is necessary to 
take into account that the values of the parameter G are 
limited by condition (|31|) : G < x 2 i/ u, rB cos 2 j3. Then, 
from Eq. (|6"2")l one can obtain the estimate (at x ~ 0.1) 
of the maximum temperatures maxTg: 



max To 



0.017 

ujr cos 2 /3 ' 



(63) 



By smooth and slow variation of the saturation param- 
eters x an d G (proportional to the intensities of the fluc- 
tuating and coherent field components) it is possible to 
implement a continuous transition from "Super-Doppler 
spatially localized states" of the atom ensemble (with the 
temperatures T ~ T sup » To) to "Sub-Doppler states" 
(with the temperatures To ~ T sub <C Td)) or to any "in- 
termediate" state with the necessary temperature value 
from the interval [T sub T sup \. T sub < T < T sup . 

However, the changes of the parameters G and X are to 
be matched (for example, by means of setting a certain 
relation between them, G = G(x)), in order not to violate 
the conditions of the stable spatial atom localization (|55|) 
and pip , at any values of the parameters mentioned. Let 
us give an example demonstrating the possibility of the 
continuous monotonic temperature change of the trapped 
atoms from T sup up to T su b without violating the condi- 
tions of the stable localization. 

Let Td = 10~ 3 K (7 w 1.31-10 8 s" 1 ), the wave number 
k ss 10 5 cm^ 1 , m w 200 amu, L\ = 1.2 cm, a\ w 0.1, 
cos j3 = 0.1, b — 1.5, and the matching law is 



G = ViS ( exp 



- 1 



(64) 



where 5 = 0.015. Thus, at X ^> S the parameter G > 1, 
and at X ^ the parameter G ~ s/3 X (which cor- 
responds to the optimal relation between G and X t° r 
reaching Sub-Doppler temperatures). 

Given in Fig. 2 are the dependencies (in a logarithmic 
scale) of the three main characteristics for the localized 
states of the atoms on the governing parameter x : the 
temperature, ln(To/To), capture velocity, ln(mw 2 /TD) 
and localization parameter, lnT (See Eq. (48), comments 
to inequality (31) and Eq. (55)). One can see that when 
decreasing the governing parameter X from the value X ~ 
0.07 (G w 2.57) to the value X ~ 0.001 the temperature 
is monotonely decreased from the value T sup w IOOTd ~ 
0.1 K to the value T sub - 0.1T D » 100 fiK (i.e. is 
decreased approximately by 1000 times, "running over" 
all the intermediate values between T sup and T sub ). The 
conditions of stable localizations (|55|) and pT|) are well 
satisfied. This, in particular, is revealed in the decay 
time of the localized states r (see Eq. (|57p) in the given 
example exceeding the characteristic diffusion time Td ~ 
1 s not less than by two orders of magnitude for any value 
X from the indicated range: r > 300^. Attention should 
also be paid to a rather high temperature value T c ~ 1 
K, corresponding to the value of the capture velocity 
v c on the right boundary of the X change interval. This 
may have important practical significance for the effective 
solution of the problem of initial atom loading into the 
optical trap under consideration. 



VI. CONCLUSION 



which can be reached in this regime of the optical field 
influence on the tripod-type atoms. 



We have developed a 3D kinetic model, describing the 
mechanical action of non-monochromatic optical fields 
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0.07 



FIG. 2. The dependence of the characteristics of the localized 
states of the atom ensemble on the governing parameter \ 
at continuous transition from "Super-Doppler states" (To 2S> 
Td) to "Sub-Doppler states" (To <C Td) and according to 
the matching law (|64[l : curve 1 — ln(Tb(x) /Td); curve 2 — 
0.2 x lnT(x); curve 3 — \n(mv^(x)/T D ); T D is the Doppler 
cooling limit, 7 ~ 1.31 -10 8 s _1 , k w 10 5 cm -1 , m = 200 amu, 
Li = 1.2 cm, a? « 0.1, cos/3 = 0.1, b = 1.5, T (xa) = To, 
Xa ~ 0.01. Inset — the details of the same curves in the Sub- 
Doppler temperature region {% ^ Xa)\ Tq(xb) = mti^Xs), 
Xs ~ 10" 4 . 



on tripod-type atoms. A distinctive feature of the pro- 
posed model is that light induced 3D atom kinetics is 
completely determined by the rectified radiation forces. 
Moreover, (under the conditions considered) the light in- 
duced atom motions in the three mutually orthogonal 
directions prove to be independent, which allows one to 
reduce the analysis of the 3D problem to the study of 
much more simple ID problems and to obtain approxi- 
mate analytical solutions of the kinetic equation for the 
Wigner distribution function of atoms /(r, v, t). 

Another important peculiarity of the considered 
scheme of the mechanical action of light on the tripod- 
type atoms is the presence of several independent govern- 
ing parameters. In particular, these are the parameters G 
and x, proportional to the intensities of the coherent and 
partially coherent optical field components, phase shifts 
£1 and 771, geometrical parameters (angular detunings Si, 
<5- and the angle (3 (see Fig. 1)) determining the optical 
field configuration. A purposeful choice of their values 
or relations between them makes it possible to set con- 
siderably different (desirable) regimes of the mechanical 
action of light on the atom motion. 

Note that, selection of the required controlling param- 
eter combination is getting significantly easier due to 
obtaining of explicit analytical expressions in the pro- 
posed model (for the RRFs and the distribution func- 
tion), which are convenient for analysis performance. 



Particularly, the developed theoretical model has been 
applied to analyze the problem of 3D optical cooling and 
trapping of the tripod-type atoms. We have made a de- 
scription of the long-lived spatially localized atom ensem- 
bles (trapped into deep light-induced potential wells) and 
demonstrated the possibility of obtaining both "Super- 
Doppler stable localized states" (with the effective tem- 
peratures T sup considerably exceeding the Doppler cool- 
ing limit Td (T sup ^> To) and high values of the capture 
velocity v c and "Sub-Doppler stable states" (with the 
temperatures T su b considerably lower than the Doppler 
cooling limit (T s „& -C Td)). We have also shown the 
possibility of continuous transition between these states 
or reaching any other intermediate state with the effec- 
tive temperature T su f, < Tjj < T sup without violating the 
localization stability. 

Finally, it is worth noting that the tripod-type quan- 
tum transitions (which have already been used in the 
experiments on laser cooling and trapping) are character- 
istic, in particular, of the atoms Na |31| (and the similar 
ones), as well as the ions 171 Yb + and 199 Hg + [32|, Un- 
interestingly, for the transitions mentioned (in iVa-like 
atoms) the usually accepted MOT model does not work; 
thus, in order to explain its work one should take into 
account the effects connected with the existence of ad- 
ditional energy levels: magnetically induced level-mixing 
effects |3~i| . One can see that the considered scheme of 
all-optical cooling and trapping of the atoms does not 
require similar effects to be taken into account. 

Regarding the above-mentioned tripod-type ions at- 
tention is to be paid to the fact that recently there 
has been some interest in the development of all-optical 
methods of ion trapping [34]. Besides, such methods can 
be used for the solution of a new interesting problem, that 
of obtaining ultra-cold electron-ion plasma with resonant 
ions and its long-term confinement [H, [35[ . 

Certainly, the application of the developed theoreti- 
cal model to certain real atoms and ions may require 
additional analysis which would take into account, for 
example (if necessary) additional repumping laser fields. 



Appendix: Explicit expressions for the terms of 
equations ([6]) and (O, describing the recoil effect 

Let us introduce an auxiliary function Qu(r, v, i) and 
consider the relations between the function p(r, v, t) and 
/(r, v, t) qu{r, v, t) (where i = x,y,z): 



Qu = / + <?« - X! 9 "' 



(A.l) 



(A.2) 



Then, the explicit expressions for the terms Aj, Bij in 
Eqs. ([5]), (O can be written in the following form: 
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A, 



hi 

2m 



<9v 



2 ov 



(A.3) 



" 2m 



<9v 



il 



2 y^ 



(A.4) 



Si 



2m 



dq 



<9v 



A, 



l^il 2 
A, 



<9v 



(A.5) 



Thus, from Eqs. ([6]), ([7]), taking into account Eqs. 
(|A.3|) - (|A.5j) . one can see that the induced dipole moment 
and the functions qij(r, v, t) (characterizing the coher- 
ence between \bi) and \a) states) are determined not only 
by the magnitude of the fields in a given point but also 
by their gradients (compare with the well-studied case of 
the two- level atom [l|); moreover, the contributions into 
this effect, connected with the coherent part of the op- 
tical field, are proportional to the gradients of the light 



I 

(Stark) shifts of the energy levels. Note that in a particu- 
lar case when the field is polarized along one of the axes 
of the Cartesian coordinate system, for example, along 
e x : Uj = 6 XJ U, Vj = S xj V, and 7 Eqs.©-© arc 
transformed into the averaged (over high frequency oscil- 
lations) generalized Bloch equations of the kinetic theory 
of RRFs for the two-level atoms in the bihromatic field 
201. 
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